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ABSTRACT 


Many satellites sent into orbit during the last few years 
are spinning to maintain stabilization, or for other reasons. 
The purpose of this paper is to determine the spin axis, a 
unit vector s°, of the satellite using ground tracking infor- 
mation. It will be shown that range rate information is ade- 
quate for this purpose. The determination of s° is neces- 
sary to check the onboard sensors once the spacecraft is in 
orbit. In the case of a slight onboard malfunction, the spin 
axis could not be determined otherwise. 
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by 

F. O. Vonbun 
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INTRODUCTION 

Often it is desirable or necessary to determine the spin axis of a rotating satellite in space 
in order to check onboard sensors or determine the change in spin axis because of a malfunction. 
In the latter case many of the experiments onboard can still be evaluated as long as the spin axis 
is known. Determining the spin axis by the method outlined here requires the radiating antenna 
to be displaced from the spin axis center. This condition is necessary since the Doppler modula- 
tion imposed on the range rate of the rotating satellite is used for the analysis. The range rate 
of a satellite in orbit can be determined with great precision (References 1 through 4); errors in 
the order of cm/sec have been achieved with CW- ranging systems. This fact will be utilized to 
determine the spin axis (unit vector s°)* of a rotating satellite. The only assumptions made here 
are that the antenna of the satellite is displaced from its axis of symmetry, that the spacecraft is 
spinning only, and that the satellite orbit can be, or has been, determined. 


RANGE RATE MODULATION DUE TO SATELLITE SPIN 

Some of our satellites are spinning in space with the antenna displaced from its spin axis.* 
This paper shows how the range rate r is modulated by the spinning motion of the spacecraft. Fig- 
ures 1 and 2 depict the vector diagram of the antenna motion in detail as well as the angle <// = arc 
cos (r° • s°) between the satellite spin axis s° and the satellite position vector r (or r°). This 
angle, as shown later, will be used for the determination of s°. Using vector notation, the antenna 
position vector 3 can be written as 

d = d Q cos coAt + d(s° x d Q ) sin coAt , (i) 


*A unit vector will always be designated with a zero superscript; that is, x 
tSee Reference 5 for centered turnstile antenna. 


x/|x|. 
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= ARC COS ( r° • S°) 
o - SATELLITE ROTATIONAL VELOCITY 

- STARTING ANTENNA POSITION VECTOR 
7 - ANTENNA POSITION VECTOR 

Figure 1— Antenna motion vector diagram. 

Also, differentiating Equation 1 with re- 
spect to time and introducing it into Equation 3 
yields, after some manipulation, 

Ar - -oxl(d 0 * r°) sin <^At + cjd(d Q * (7 xs°))cosoAt 

+ ( 4 ) ( 3 ° ' 7 ) “ (t) ’"( 3 ° • 7°) . (4) 

Inspecting Equation 4 shows that as long as 

« co and ~ « go , (5) 

the third and fourth terms of Equation 4 can be 
neglected.* 


where 3 ° = d • 3 ° is the starting position of the 
spacecraft antenna at t = t 0 or At = 0. Equa- 
tion 1 follows directly from Figure 1, as can be 
seen. Projecting 3 onto the range vector 7 re- 
sults in: 

Ar = (d • r°) = 7 - (d • r) . (2) 

Since the spacecraft is spinning, the value 
Ar changes with time (and with orbital motion, 
which is assumed to be negligible as shown 
later); the following is obtained by differentiat- 
ing Equation 2 with respect to time: 

Ar = [(d • r) + (d ■ r)] - (3 • 7) rj • (3) 


/ 



/ 
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RANGE RATE r 
IS MODULATED BY o, 

THE ROTATION OF THE SATELLITE 

Figure 2— Satellite and antenna details where d equals 
antenna distance from spin axis. 


♦Example: jr*f = |v | (neglect earth rotation) - 8000 m/sec 

n = 9 revolutions per min.; co - 27rn/60 radians/sec = 1 sec -1 (approx.) 
r - 400 km (slant range) 
than: |r*|/r - 2.10" 2 sec -1 « 1 sec -1 . 
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Since 3° is an arbitrary starting vector 
(unit vector) it will be chosen so that 

(d 0 ° -7°) = 0 , (6) 

or 

do = 7IH^r(r° x s°) . 

Introducing Equation 6 into Equation 4, neglect- 
ing the terms just mentioned, yields 

Ar = (cod sin t//) cos o;At , (7) 

(with v/r « go, r/r«co) as the range rate mod- 
ulation caused by the spinning motion of a 
satellite. 

If the range rate modulation Ar. from the 
i th ground tracking station is given as shown 
in Figure 3, the angle > p k between the spin axis 
and the position vector r. (or r ° ) can be de- 
termined from Equation 7; that is, 



since co = 2n/T is given by the same process. A 
better determination of cj can be obtained from 
the power spectrum of the range- rate informa- 
tion as shown in Figure 4. Angle </>. will now be 
used to determine the components of the unit 
vector s°, the satellite spin axis. 

DETERMINING COMPONENTS OF SATELLITE 
SPIN AXIS UNIT VECTOR s° AND 
THEIR ERRORS 

Adapting the legend in Figure 1 (see also 
Figure 5) gives: 

cos \Jj . = (r ° • s°) , (9) 



WL I w in IYIVIIVM 

(a) " Spin modulated" range rate TIME 



( b) Range rate modulation ( orbital motion subtracted) 


Figure 3— Range rate modulation. 
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Figure 5— Satellite spin axis in relation to ground- 
station line of sight. 


where r? is the unit satellite position vector of 
the i th ground station. From Equation 8, the 
angle can be determined; thus only s° is un- 
known in Equation 9. A straight-forward solu- 
tion of Equation 9 would be to measure the Ar. 
from three stations and determine the compon- 
ents s 13 s 2 , s 3 from Equation 9. 

However, a much better solution can be ob- 
tained by using all available information; that is, 
by using much more than three measurements. 

In practice, many ground stations take a large 
number of range rate measurements. This in 
turn dictates a least- square solution, since 
more than the necessary three equations (meas- 
urements) are available. 

The dot product of Equation 9 can be writ- 
ten also in component form as 

cos./;. = i-j, s, + r ;2 s 2 + r i3 s 3 , (10) 


where the values r. l9 r i2 , r i3 and s 13 s 2 , s 3 are the components of the unit vectors and s° 3 
respectively. 

Using many measurements of Ar. and a> i from which 0. can be calculated from Equation 8, one 
has a set of linear equations at hand. From Equation 10 one obtains for n measurements of Ar. 
and the following: 


a i - cos \p i = r , , s, + r 12 s 2 + r 13 s 3 
a 2 = cos 0 2 
a = cos 0 

n 

Equation 11 suggests a matrix treatment which will be followed from this point on. In matrix 
form Equation 11 reads 


r 2 1 S 1 + r 22 S 


2 + r 23 


S 1 + r n2 S 2 + r n3 S 3 


(ID 


nX 1 ) nX 3 ) ' ^(3X1) 


( 12 ) 


where the matrices A, R, and s (the unit vector of the satellite spin axis) are as follows: 



a i 


r ll r 1 2 r 1 3 


V 

A = 

a 2 

, R = 

F 2 1 r 2 2 r 2 3 

, s = 

S 2 



( Axl) 

r n2 r n£ 

( nX3) 

S 3_ 


mil 


I IBI 


III 



To solve Equation 12 in the least- square sense, standard procedures are being applied (Ref- 
erences 6 through 12). Since in practice Equation 12 as such will never exist, because of en- 
hanced errors in the measurements as well as in the analyses, one must write instead: 


A(nXl) nX3) 


3X 1 ) + € ( nX 1 ) 


( 13 ) 


where e is an error matrix with expectation 


E(£) = 0 


(14) 


and the dispersion matrix 


V(s) 


( nX n ) 


E(- T ) 


( nX n ) 


(15) 


where i (nXn) is a unit matrix and cr a 2 . is the variance of the measurements made. This means 
that in each of Equations 11 an error term must be added to correspond with actual conditions. 

Equation 13 must be used to determine the matrix S or the unit vector s°. The principle of 
least squares requires that S be solved for in such a way that the sum of the errors squared (e T e) 
is a minimum (e T is the transpose of e). A necessary condition for this is that 

(e T e) = 0 ; (16) 


or, using Equation 13, one obtains from Equation 16 

Js [(a t -S t R t ) • (A-RS)] = 0, 

which yields, after some calculations, the best estimate for S; that is, 

^(3Xl) ~ (^ T ®)( 3X 3 ) 3Xn) \ nX 1 ) * 


(17) 


Equation 17 represents the solution to the problem of determining satellite spin axis. 

In order to estimate the errors made in the determination of s°, its dispersion matrix must 
be calculated. From Equations 17 and 15 one obtains: 


S = (r t r)” 1 R t (RS + e) . 


(18) 
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Using Equation 15 and introducing Equation 18, one obtains, 


V(S) = E(€£ t ) = erf. (r t r)~ 1 


( 19 ) 


for the dispersion matrix of S. 

ESTIMATING MEASUREMENT VARIANCE 

In order to evaluate the dispersion matrix V(S) (3X1) the measurement variance must be esti- 
mated also. From Equation 8 Ai-j and gj. are measured quantities which in turn have errors as- 
sociated with them. These can be estimated from the noise of the data (this is to be assumed here 
as known). From Equations 8 and 11 one can write 

1 /^imax\ 

sin ^i - d \ o>. J* (20) 


and 


cos \b . - a. . 

r X 1 

Varying both equations, one obtains for the variation in a. 


Sa . 



(21) 


Assuming now that the errors are uncorrelated and replacing the value Ai-. max by (Ar. max + crAr imax ) — 
since it can never be zero, because of the inherent errors in its determination— one obtains from 
Equation 21: 


Cr 


a i 



( 22 ) 


With the measurement variance o- a \ and the matrix R known, the disposition matrix V(S) given by 
Equation 19 can be evaluated. Equation 22 shows also that the limitations for cr ai ; (Ar imax 
+ crAr imax /Wi d) must be smaller than 1 to prevent excessive errors. 
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CONCLUSIONS 


It has been shown how the range- rate modulation of a spinning satellite with a displaced an- 
tenna can be used to determine its spin axis (unit vector) in space. An estimate of the spin axis 
errors associated with this method has also been presented. It is interesting to note that no addi- 
tional equipment is necessary on the satellite or on the ground. This method may well be suited 
for determining the spin axis of communications and navigation satellites. If necessary, it could 
also be used to check the axis of spin- stabilized kick stages by adding a Doppler transponder and 
antenna to the stage. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, July 27, 1967 
550-311-07-21-51 
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